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ABSTRACT 

We present Spitzer/IRS mid-infrared spectra for 15 gravitationally lensed, 24 /^m-selected galaxies, 
and combine the results with 4 additional very faint galaxies with IRS spectra in the literature. 
The median intrinsic 24 /zm flux density of the sample is 130 /iJy, enabling a systematic survey of 
the spectral properties of the very faint 24 /im sources that dominate the number counts of Spitzer 
cosmological surveys. Six of the 19 galaxy spectra (32%) show the strong mid-IR continuua expected 
of AGN; X-ray detections confirm the presence of AGN in three of these cases, and reveal AGNs in two 
other galaxies. These results suggest that nuclear accretion may contribute more flux to faint 24 /im- 
selected samples than previously assumed. Almost all the spectra show some aromatic (PAH) emission 
features; the measured aromatic flux ratios do not show evolution from z = 0. In particular, the high 
S/N mid-IR spectrum of SMM J163554. 2+661225 agrees remarkably well with low-redshift, lower- 
luminosity templates. We compare the rest-frame 8 /xm and total infrared luminosities of star-forming 
galaxies, and find that the behavior of this ratio with total IR luminosity has evolved modestly from 
z=2 to z=0. Since the high aromatic-to-continuum flux ratios in these galaxies rule out a dominant 
contribution by AGN, this finding implies systematic evolution in the structure and/or metallicity of 
infrared sources with redshift. It also has implications for the estimates of star forming rates inferred 
from 24 /im measurements, in the sense that at z ~ 2, a given observed frame 24 /im luminosity 
corresponds to a lower bolometric luminosity than would be inferred from low-redshift templates of 
similar luminosity at the corresponding rest wavelength. 
Subject headings: galaxies — infrared: galaxies: high-redshift — galaxies: evolution 
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L INTRODUCTION 

The Spitzer Space Telescope (Werner et al. 2004) has 
been tremendously successful at detecting star-forming 
galaxies and active galactic nuclei (AGN) at z > 1 by 
their emission in the observed-frame 24 /Ltm band. Down 
to the confusion limit of ~ 50 /iJy at 24 /j,m (Dole et 
al. 2004), the MIPS instrument (Rieke et al. 2004) can 
detect luminous infrared galaxies (LIRGs)^ out to z ~ 2, 
ultra- luminous infrared galaxies (ULIRGs)^° out to z ^ 
3, and hyper-luminous infrared galaxies (HLIRGs)^^ out 
to even higher redshifts. 

Detectability at z > 0.7 is increased by strong aro- 
matic emission features, at rest-frame 6-12 /im, that pass 
through the 24 /im band. Indeed, photometric redshifts 
place ~ 25-30% of the faint 24 /im sources at redshifts 
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above 1.4 (Le Floc'h et al. 2005; Perez-Gonzalez et al. 
2005; Caputi et al. 2006; Wang 2006). 

For extremely high-luminosity galaxies, spectra have 
been obtained with the Infrared Spectrograph (IRS, 
Houck et al. 2004), e.g. Houck et al. (2005); Yan et 
al. (2007); Menendez-Delmestre et al. (2007). However, 
since spectral properties depend strongly on luminosity, 
results obtained for hyper- luminous galaxies may have 
limited applicability for the bulk of the IR-detected pop- 
ulation. Spectra of lower-luminosity galaxies have been 
obtained only in a few cases (e.g. Lutz et al. 2005; Teplitz 
et al. 2007). 

Numerous authors (e.g. Caputi et al. 2006; Choi et al. 
2006; Reddy et al. 2006) have used the observed 24 /im 
band as an estimator of star formation rates. Since most 
of the infrared power is radiated redward of the observed 
band, the 24 /zm diagnostic must be calibrated using 
low-redshift templates (e.g. Chary & Elbaz 2001; Dale 
& Helou 2002; Lagache et al. 2004; Brandl et al. 2006; 
Armus et al. 2007). The central assumption in these 
works is that the spectra and spectral energy distribu- 
tions of high-z galaxies are matched accurately by the 
low-redshift templates. For 0.7 < z < 3, the 24 /im band 
fluxes are strongly influenced by aromatic band emission. 
Evolution over time either in the behavior of these aro- 
matic bands, in metallicity, or in the geometry and ra- 
diative transfer within these infrared sources, could all 
undermine this central assumption. 

To probe these possibilities, we have been obtaining 
IRS spectra for intrinsically faint 24 /im-selected galax- 
ies at 1 < z < 3. Our targets are strongly lensed by the 
gravitational potential of foreground clusters of galax- 
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ies, such that their fluxes are amphfied by factors of 3- 
25. Though these galaxies have observed 24 /im flux 
densities of '--^ 1 mJy, sufficient for useful, efficient spec- 
troscopy with IRS, their intrinsic luminosities are down 
in the LIRG to ULIRG range. 

In this paper, we present IRS spectra for the first 15 
galaxies in the sample. The median observed flux den- 
sity is fi/(24 /xm) = 0.82 mJy, but the median intrinsic 
flux density is nearly an order of magnitude lower. We 
compare their aromatic feature flux ratios, and the ra- 
tios of their X-ray, aromatic, and total infrared luminosi- 
ties to those of low-redshift LIRGs and ULIRG samples. 
We examine in detail the high-quality spectrum of SMM 
J163554.2+661225, a LIRG at z = 2.516. We merge our 
sample with four additional objects from the literature 
(Lutz et al. 2005; Teplitz et al. 2007) to provide a first 
systematic look at the IRS-measured properties of the 
faint 24 ^m galaxy population. 

We assume ^m = 0.265, flA = 0.735, and Hq = 
72 km s~^ Mpc"'^ (Spergel et al. 2003; Freedman et al. 
2001). 

2. SAMPLE SELECTION 

Our goal was to select sources with intrinsic flux den- 
sities close to, or even below, the MIPS 24 /im confusion 
limit. We first obtained deep 24 iim images of galaxy 
clusters, generally identified as being massive through 
high X-ray luminosity. From these 24 /im images, we 
selected sources within 1.5' of the cluster center that: 

• had observed 24 /im fiux densities above 0.4 mJy 
(to make efficient use of IRS). 

• had faint, irregular, or arc-like HST counterparts. 
This excluded galaxies at the cluster redshift.^^ 

• had known spectroscopic or probable photometric 
redshifts > 1.0. 

• had lensing amplification factors calculated to be 
greater than 3. This restricted selection to clus- 
ters with good lensing models. Since amplifica- 
tion depends on the source redshift, if the red- 
shift was unknown, we calculated the amplification 
for 1 < z < 3 and required amplifications factors 
greater than 3 over that range. 

Since submillimeter sources in lensing clusters have 
been studied extensively, with published amplifications 
and redshifts, we prioritized sources in the target list 
with published sub- mm detections. Our sample contains 
15 galaxies, targeted in GTO campaigns 82 and 30775. 
Eight of these galaxies have sub-mm detections in the lit- 
erature, as referenced in Table 1. Figure 1 shows optical- 
band images of the sample galaxies. 

We supplement the sample with four sources from the 
literature. We use three of the four sources from Teplitz 
et al. (2007), with /^(24 /im)= 0.13-0.2 mJy; these are 
unlensed galaxies at 1 < z < 3 with very long (9-12 hr) 
IRS exposures, and are not known sub-mm sources. We 

^^ A2219a and A2261a lack HST imagery; we instead used 
ground-based F606W images from the Steward Observatory 90 inch 
telescope. Both galaxies are known sub-mm sources (see Table 1), 
and can thus be added to the sample despite the lack of high- 
resolution optical images. 



also use the lensed z — 2.81 sub-mm source in A3 70 from 
Lutz et al. (2005)" with f^{24 /<m)=1.36 mJy 

Source coordinates, redshifts, amplifications, and flux 
densities or our sample and for the four literature sources 
are all listed in Tables 1 and Table 3. 

3. OBSERVATIONS AND DATA REDUCTION 

This paper is based on a combination of imaging with 
MIPS, spectroscopy with IRS, archival imaging from 
Chandra, and submillimeter photometry from the liter- 
ature. Exposure times are summarized in Table 1. 

3.1. Spitzer/IRS spectroscopy 

Low-resolution IRS spectra were obtained as part of 
Spitzer GTO programs 82 and 30775 (PI G. Rieke). 
Data reduction used a package developed by D. Fadda, 
which has already been successfully used to reduce low- 
resolution IRS spectra of faint high-z sources (Yan et 
al. 2007). Residual background, rogue pixels^** (pixels 
with dark current values abnormally high and variable 
with time) and cosmic ray hits have to be corrected to 
obtain an optimal reduction of the data. The package 
considers all the bidimensional frames produced by the 
IRS/SSC pipehne. First, background and noise images 
are produced by coadding frames after masking target 
and serendipitous spectra on each frame. (The back- 
ground is better estimated for the PID 30775 observa- 
tions where both LLl and LL2 orders provide redundant 
background measurements, compared to the PID 82 ob- 
servations using only LLl, where the background can 
only be determined by differencing nods and thus suffers 
from contamination by faint sources.) 

A robust statistical estimator (biweight) is used to min- 
imize the effects of deviant pixels on the coadded value. 
This part is iterative in order to allow manual identifi- 
cation of sources. Rogue pixels are then identified by 
computing the dispersion of the noise around every pixel 
and flagging pixels which are > 5a deviant from the mean 
local value. Spectral extraction is done optimally (e.g., 
by using the PSF proflle to weight the spectrum), taking 
into account the spectral distortion, and also rejecting 
pixels affected by cosmic rays. Finally, the spectra ob- 
tained at each position are co-added. 

We downloaded and reduced one archival IRS spec- 
trum, from program 3453 (PI P. van der Werf). This 
is SMM J163555. 2+661238, image "A" in Kneib et al. 
(2004), one of three images of the same triply-imaged 
sub-mm galaxy in Abell 2218. We obtained an IRS 
spectrum of the brightest image of this galaxy, SMM 
J163554.2-h661225 (image "B" in Kneib et al. 2004); 
its 24 /xm flux is 1.8 times brighter than the SMM 
J163555. 24-661238. In the archival spectrum, the source 
was mis-centcrcd across the short axis of the slit, which 
reduced throughput. Because its signal-to-noise ratio is 
low, we use the archival spectrum only to confirm the 
detection of features, and use the spectrum from our pro- 
gram for all analysis. 

The IRS spectra are presented in Figure 2. Spectra 
obtained in programs 82 and 30775 for the same object 

^■^ The two other sources from those papers, "2-x" from Teplitz 
et al. (2007) and the Abell 2125 source from Lutz et al. (2005), 
have insufficient signal-to-noise for our analysis. 

^"^ http://ssc.spitzer.caltech.edu/irs/roguepixels/ 
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are plotted separately; observations from the latter pro- 
gram have superior background subtraction. Flux den- 
sities are as observed, without any correction for gravi- 
tational amplification. When redshifts were known from 
the literature, or are evident from aromatic emission fea- 
tures in the IRS spectrum, then the upper x-axis shows 
rest wavelength. Only one source (A2261a) lacks a red- 
shift. Redshifts are listed in Table 1, with a reference if 
the redshift is from the literature. 

The rest-frame spectra were run through PAHFIT 
(Smith et al. 2007) to fit the continuum and aro- 
matic (PAH) features simultaneously. Measured aro- 
matic fluxes are reported in Table 1. 

We allowed PAHFIT to fit the silicate absorption (e.g., 
we did not specify "NO.EXTINCTION"). Fixing the sil- 
icate optical depth at zero changes the best-fit 7.7 /im 
feature fluxes by 4% (median), the best-fit 8.6 /um fea- 
ture fluxes by 8% (median), and the best-fit 11.3 /im 
feature fluxes by 11% (median). 

We report fluxes in table 1 for features in which the 
measured flux is at least twice the uncertainty calculated 
for that line. We add back in a few features, mostly at 
11.3 /xm, which are obviously detected even though the 
errorbars are large. 

3.2. Spitzer/MIPS 24 fJ-m and 70 jim photometry 

MIPS photometry-mode images at 24 /im and 70 /im 
were obtained through Spitzer GTO program 83 (PI 
G. Rieke). Exposure times are given in Table 1 and pho- 
tometry is given in Table 3. The data were reduced and 
mosaicked using the Data Analysis Tool (Gordon et al. 
2006) with a few additional processing steps (Egami et al. 
2006). 

Photometry at 24 /tm was obtained by PSF fitting, 
using the IRAF implementation of the DAOPHOT task 
allstar. The PSF was created empirically using all avail- 
able images from program 83. An aperture correction of 
1.131 was applied, calculated from a Tiny Tim model of 
the 24 /im Spitzer PSF that extends to r = 220" (C. En- 
gelbracht, priv. comm.)^^ 

Aperture photometry at 70 /im was performed on the 
1.0 pixel scale images. The aperture radius was 16" 
(twice the half- width at half-max), the background an- 
nulus had radii of 18" and 39" (the position of the first 
airy ring), and an aperture correction of 1.968 was ap- 
pfied.^^ If the source was undetected, we took the upper 
limit as 2x the la sky noise plus any positive fiux in the 
source aperture. 

3.3. Spitzer/IRAC photometry 

IRAC images at 3.6, 4.5, 5.8, and 8 /im were obtained 
as part of GTO program 83 (PI G. Rieke) for five of 
the six clusters; IRAC data for AbcU 1835 were obtained 
through GTO program 64 (PI G. Fazio). Exposure times 
are listed in Table 1. IRAC images were mosaicked as 
described by Huang et al. (2004). 

Photometry was obtained by aperture photometry 
with sky subtraction. The high density of cluster galax- 
ies in the cluster cores necessitated use of irregularly- 
shaped polygons for both target and sky apertures. The 

^^ The aperture correction is specific to the DAOPHOT param- 
eters used, in this case sky annulus radii of 31.3" and 40", and a 
PSF defined out to r = 22.41". 

^^' http://ssc.spitzcr.caltcch.edu/rnips/apcrcorr/ 



flux density within these polygonal regions was deter- 
mined using CIAO^^ Versions 3.2.1 and 3.4.1.1. IRAC 
photometry is listed in Table 3; the values do not include 
aperture corrections. 

3.4. HST photometry from WFPC2 and ACS 

From the HST archive, we downloaded all publi- 
cally available ACS images for the clusters. We used 
Multidrizzle (Koekemoer et al. 2002) to distortion- 
correct each flat-fielded pipeline image, and then cross- 
correlated the images to measure the small (< 1") coor- 
dinate registration offsets with high precision. The mea- 
sured offsets were then used by multidrizzle for the final 
mosaicking. The data were photometcrcd using the same 
technique as for the IRAC bands. 

We also used reduced WFPC2 images kindly made 
available by D. Sand, which are described in Sand et 
al. (2005). 'ACS and WFPC2 cutouts of the galaxies in 
our sample are shown in Figure 1. 

3.5. Chandra ACIS imagery 

From the Chandra archive, we downloaded all publi- 
cally available ACIS observations of clusters with IRS 
targets in this sample: observation ID numbers 529 and 
902 for MS0451; 1663, 5004, and 540 for A1689; 495 and 
496 for A1835; 1454, 1666, and 553 for A2218; 896 for 
A2219; 5007 and 550 for A2261; 4193, 500, and 501, for 
A2390; and 1562 for AC114. 

We used CIAO (3.3) to update the CTI correction 
when necessary and remove data obtained during periods 
of high solar flares. We mosaicked multiple overlapping 
observations following the CIAO thread "Reprojecting 
Images: Making an Exposure-corrected Mosaic" . 

For X-ray detections, fluxes were determined from the 
individual (not mosaicked) observations as described in 
Marcillac et al. (2007). X-ray flux upper limits (3(t) 
were determined from the mosaicked images and mo- 
saicked exposure maps as in Donley et al. (2005) , except 
that the 90% encircled energy radius was used, and the 
Monte-Carlo sky apertures were taken close to the source 
(within 30") in order to sample a representative cluster 
background. 

Chandra X-ray fluxes and upper limits for the 2-8 keV 
and 0.5-8 keV bands are reported in Table 1. 

3.6. Long-wavelength photometry from the literature 

Eight galaxies in our sample, and one galaxy from the 
extended (literature) sample, have submillimeter detec- 
tions in the literature; these are listed and referenced 
in Table 1. Four sources were detected at 15 /im by 
ISO: A851a and A2218b (Barvainis et al. 1999), A2390b 
Lemononet al. (1998), and A2390c (Metcalfe et al. 2003). 
We will use this photometry to help determine the total 
infrared luminosities of the sources. 

4. DISCUSSION 

4.1. Aromatic Feature Flux Ratios in High-Redshift 
Galaxies 

We now test whether the mid-IR spectra of star- 
forming galaxies, with their prominent aromatic emission 

^"^ http://cxc.harvard.edu/ciao/ 



Rigby et al. 



features, show evidence for evolution from z = to high 
redshift. 

Given its high quahty, we first examine the IRS spec- 
trum for SMM J163554.2-t-661225 in detail The image 
is amphfied by a factor of 22 ± 2 (Kneib et al. 2004), 
such that without lensing, its 24 /im flux density would 
have been 53 /zJy, at the MIPS confusion limit. The 
galaxy redshift, from Ha, is z = 2.5165 ± 0.0015 (Kneib 
et al. 2004). Its IRS spectrum, presented in Figure 3, 
shows highly-significant detections of aromatic features 
at 6.22, 7.7, 8.33, and 8.61 /im. (The lower-quality spec- 
trum from program 3453 confirms the detection of these 
aromatic features, though the two longer-wavelength fea- 
tures are detected at low significance.) Overplotted with 
arbitrary normalization are two spectral templates: a) 
the starburst galaxy NGC 2798 from the SINGS survey 
(Dale et al. 2006); and b) the average spectrum of 13 
nearby starburst galaxies with little apparent AGN con- 
tribution (Brandl et al. 2006), whose mean IR luminosity 
is 4.9 X 10^° Lq. Coincidentally, this is also the L(TIR) 
of NGC 2798. 

The z — 2.516 spectrum in Figure 3 closely resem- 
bles the low-redshift, lower-luminosity (by a factor of 15) 
templates. Aromatic flux ratios are reported in Table 1, 
along with those we measure for the Brandl template, 
and those reported for the SINGS sample. The strength 
of the 8.6 ^m feature in SMM J163554.2+661225 rela- 
tive to those at 6.2 /im and 7.7 /im, is within the 10%- 
90% variation within the SINGS sample. In addition, the 
flux ratios of SMM J163554.2-h661225 are a close match 
to those of the Brandl et al. (2006) average starburst 
template; while 7.7 /im is relatively stronger than the 
template, the difference is comparable to the standard 
deviation of line ratios for individual galaxies that were 
averaged to create the template. Thus, we conclude 
that this z = 2.5 galaxy has aromatic feature flux 
ratios that are consistent with those observed for 
lower luminosity, z = starbursting galsLxies. 

While the other spectra in the sample have lower 
signal-to-noise ratios, they are still sufficient, in the 
aggregate, to test whether high-redshift galaxies have 
markedly different aromatic flux ratios. Figure 4 com- 
pares, for four aromatic feature flux ratios, the observed 
ratios in our sample with those measured for the low- 
redshift Brandl et al. (2006) template. Though the er- 
rorbars in individual measurements are sizable. Figure 4 
confirms the result seen for SMM J163554.2+661225: 
star-forming galaxies at 1 < z < 3 have aromatic fea- 
ture ratios consistent with those of z = star-forming 
galaxies. 

Variations in aromatic feature ratios have been re- 
ported in certain HII regions — this may be attributable 
to variations in size, composition, and ionization state of 
the carriers (Draine & Li 2007 and references therein). 
In particular, large PAHs emit more strongly at 11.3 /xm, 
and neutral PAHs emit much more strongly at 3.3 and 
11.3 /im than charged PAHs (Draine & Li 2007; Galliano 
2006; AUamandola et al. 1989). That our measured fea- 
ture ratios at z > 1 are consistent with those measured 
at z ~ suggests that the dust grain size and ionization 
distributions are not strongly evolving. 

4.2. Evidence for Compact Source Accretion 



Active galactic nuclei are detectable in multiple ways — 
by the dust they heat that radiates in the mid-IR; by 
X-rays emitted from the accretion disk coronae; and by 
high-excitation or broad emission lines. We now examine 
the evidence for AGN activity in our extended sample. 

The first AGN diagnostic we examine is X-ray lumi- 
nosity, plotted in Figure 5 against aromatic luminos- 
ity. Most of the sample are X-ray non-detections, at 
limiting luminosities that rule out X-ray-loud QSOs or 
bright Seyferts, although an X-ray-weak or highly ob- 
scured AGN could still be present (c.f. Donley et al. 
2005; Alonso-Herrero et al. 2006). Two of the four 
galaxies from the literature, and three of the galaxies 
in our program are detected in X-rays: two have no ap- 
parent aromatic features (sources A2261a and A2390a); 
and the other, source A2390b, is interesting in that it 
contains a luminous X-ray-emitting AGN, yet its spec- 
trum still shows aromatic features. That source emits 
roughly equal power in the 7.7 /im aromatic feature 
(scaled from the 11 /im feature) and at 10-30 keV, remi- 
niscent of the more luminous z — 1.15 source CXO GWS 
J141741. 9-^522823 discussed by Le Floc'h et al. (2007). 

We now consider a second AGN diagnostic, the rela- 
tive contribution of aromatic versus continuum emission 
to the mid-IR fiux. At low redshift, low aromatic feature 
equivalent widths have been demonstrated as an effective 
AGN diagnostic (Armus et al. 2007; Brandl et al. 2006, 
and Devost et al. (in prep.)). However, it is extremely 
difficult to measure equivalent widths accurately for the 
high-redshift IRS sources, due to the limited wavelength 
baseline, the fact that the aromatic features have broad 
wings, and most importantly, imperfect sky subtraction 
which adds a (positive or negative) pedestal to each spec- 
trum. For example, for SMM J163554.2-h661225 (our 
highest-quality spectrum), we cannot measure an accu- 
rate equivalent width due to the difficulty in determining 
the true continuum level. We are able to set a lower limit 
of > 0.4 /im (rest-frame, conservatively assuming a high 
continuum level) , which is a typical value for low-redshift 
star-forming galaxies. For the other spectra in the sam- 
ple, equivalent widths are even harder to measure. 

Therefore, we instead create an aromatic-to-mid-IR 
flux ratio, which we deflne as the ratio of the flux in the 
7.7 /im PAH feature (in erg s"^ cm"^^ ^s fit by PAHFIT) 
to the MIPS photometry at 24 /im (in Jy, corrected for 
bandwidth compression to z = 0.) This metric does not 
suffer as strongly from the difficulty in determining the 
continuum level. In the low-redshift, star-forming com- 
parison sample we construct in § 4.4, this ratio ranges 
from 7-32 xl0~". We thus take 7 x 10"" as a dividing 
line between spectra dominated by star formation (above 
the value), and spectra with substantial AGN contribu- 
tion (below the value). This discriminant selects the fol- 
lowing mid-IR spectra^* as having substantial AGN con- 
tribution: A370a, MS0451a, A1689a, A2218c, A2261a, 
A2390a. Thus, of 19 galaxies with adequate spectra, 
6 show strong indications of AGN contribution to their 
mid-IR outputs. 

A third indication of AGN activity is a mid-IR spec- 
trum that rises steeply with increasing wavelength. Four 

^* For A2390b, we scale from the flux in the 11 /im aromatic 
feature, using the average 7.7/11 /im flux ratio of 3.6 from Smith 
et al. (2007). 
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spectra (A2261a, A2218c, A1689a, and A370a) show this 
behavior; none has detected aromatic emission. 

A fourth indication of AGN activity is the presence 
of broad or highly excited lines in an optical spectrum. 
Of the spectroscopy that has been published, or is in 
preparation, for our sample (see references in Table 1), 
three show evidence for AGN: A1689a has a Keck spec- 
trum showing highly ionized neon (J. Richard et al. in 
prep.); A2390a has a Keck spectrum (see Appendix) 
which shows Lyman a with FWHM ^ 945 km s~^, typ- 
ical of an AGN narrow line region; and the spectrum 
of A370a is reported to contain AGN lines (Lutz et al. 
2005). 

A fifth indication of AGN activity is a high ratio of 
[Ne III] 15.5 /xmto [Ne II] 12.8 /xm. Solar-metaUicity 
starbursts are not observed to excite this line ratio 
above unity (Thornley et al. 2000; Rigby & Rieke 2004), 
whereas the harder continuum of AGN do. Eight spec- 
tra cover rest-frame 12.8/im; [Ne II] is detected in 6 of 
these. Six spectra cover rest-frame 15.5 /im; [Ne III] 
is detected in 4. Of the four galaxies where both lines 
are detected (A1689a, Ab2218b, A239Gb, and A2667a), 
only in A1689a does the [NeIII]/[NeII] ratio exceed unity: 
2.15±0.76. This high line ratio offers additional evidence 
that A1689a hosts a luminous AGN. The other, lower line 
ratios are more difficult to interpret. 

The results for these AGN diagnostics are compiled in 
Table 1. Combining these diagnostics, 8 galaxies out of 
the sample of 19 have at least one indication of AGN 
activity (X-ray detection, low aromatic contribution to 
the mid-IR flux, rising mid-IR spectrum, or optical AGN 
lines). Five galaxies have two different diagnostics that 
indicate AGN. Thus, we find that 42% of the faint 24 /zm 
sources in this sample have some nuclear activity; the 
95% confidence interval is 23%-63%. Further, in 32% of 
the sample, the AGN strongly affects or dominates the 
mid-IR spectrum; the 95% confidence interval is 15%~ 
53%. (Confidence intervals were calculated using eqn. 26 
of Gehrels (1986).) 

The high proportion of AGN in our sample of high- 
redshift LIRGS and ULIRGs mirrors the general behav- 
ior of local IR-luminous galaxies, where AGN become 
substantially more common with increasing luminosity 
(e.g., Lutz et al. 1998). Though the errorbars are large 
given the small sample size, the high proportion (32%) of 
AGN sufficiently powerful to dominate the mid-IR spec- 
tra may be inconsistent with the frequent assumption 
that ~ 100 /iJy, z > 1 sources are dominated by star 
formation, with AGN a modest contaminant. Our re- 
sult, if confirmed with better statistics, has implications 
for 24 /im-derived star formation rates in deep surveys, 
where it is often assumed that the mid-IR outputs of 
high-redshift, faint galaxies are dominated by strong aro- 
matic features from star formation. It also implies that 
a higher fraction of the integrated 24 fim fiux on the sky 
may arise from accretion than previously thought. 



4.3. Total infrared luminosity 

We now determine the total infrared (TIR) luminosi- 
ties of the lensed sources, so that we may examine the 
dependence of aromatic feature emission on total infrared 
luminosity, and how this behavior may evolve with red- 
shift. 



4.3.1. SMM J163554.2+661225 

The sub-mm galaxy behind Abell 2218 has exception- 
ally good photometric coverage from 0.4 to 850 fim. 
Since it is triply-lensed, we increase the signal-to-noise 
ratio of its multiband photometry by using fiuxes from 
all three images, weighed by the amplifications given in 
Kneib et al. (2004). (Bright neighboring galaxies forced 
us to use only image A and B for the optical and IRAC 
photometry.) Figure 6 plots this photometry. 

We fit redshifted Chary & Elbaz (2001) 'and Dale & 
Helou (2002) templates to the 24, 70, 450, and 850 iim 
photometry, corrected for lensing amplification. For 
each template family, we determine the best-fit template 
and the range of acceptable templates, then determine 
L(TIR) for those templates in a method consistent with 
the definition given in Table 1 of Sanders & Mirabel 
(1996) using 12, 25, 60, and 100 y^m photometry. We 
estimate the intrinsic total infrared luminosity of SMM 
J163554.2-t-661225 as: 

• 7.6 ± 1.9 X 10" Lq for Dale & Helou (2002) tem- 
plates. 

• 1.3to:2 X 10^^ Lq for Chary & Elbaz (2001) tem- 
plates. 

As Figure 6 shows, the best-fit Chary & Elbaz (2001) 
template over-predicts all > 24 /im points, especially at 
70 /im. The best-fit Dale & Helou (2002) template is 
a better match to the shape of the overall infrared SED 
and also has high aromatic feature equivalent widths, 
consistent with the spectrum of the source. Thus, this 
galaxy's intrinsic luminosity is that of a LIRG. 

4.3.2. Other Galaxies 

We adopt a similar approach for the other galaxies 
with sufficient far-infrared and submillimeter data to 
provide reasonable constraints. We tried both Dale & 
Helou (2002) and Chary & Elbaz (2001) templates, fit- 
ted to fall just below any upper limits so long as they 
were compatible with the other constraints, i.e., to give 
a maximal luminosity. Two galaxies are fit poorly by a 
star-forming template, and much better by an AGN tem- 
plate, Mrk 231; MS0451a may be, as well. The ranges of 
acceptable luminosities and the preferred templates are 
listed in Table 3. Dale & Helou (2002) templates pro- 
vided generally better fits to our sample than Chary & 
Elbaz (2001) templates, as expected given the results of 
Marcillac et al. (2006) for < z < 1 galaxies, fitting ISO, 
Spitzer, and VLA photometry. 

4.4. 8 iim (rest) -to- Total Luminosity Ratio 

We now consider the aromatic and total infrared (TIR) 
luminosities of the star-forming galaxies in the lensed 
sample, and how this ratio behaves as a function L(TIR). 
We compare to higher-luminosity star-forming galaxies 
at z ~ 2, as well as comparable-luminosity star-forming 
galaxies at z ~ 0.8 and z ~ 0. We interpret the prop- 
erties of the high redshift galaxies strictly in terms of 
properties that are well-determined for local analogs, so 
we can test whether inconsistencies emerge. 

4.4.1. Ttie z ~ comparison sample 



6 



Rigby et al. 



To probe rest-frame 8 fim and total infrared (TIR) 
luminosities for local galaxies, we use the SINGS sam- 
ple (photometry from Dale et al. 2005) supplemented 
in the ULIRG range by IRAS 09111-1007, 10565-^2448, 
12112+0305, 14348-1447, 17208-0014, 22491-1808, and 
Arp220 (Farrah et al. 2003; Armus et al. 2007), all 
selected because they appear to be dominated by 
young stars rather than AGN.^^ We also include IRAS 
00262+4251, 01388-4618, 02364-4751, 16474+3430, 
23128-5919, and 23365+3604 from Rigopoulou et al. 
(1999). The strength of the aromatic features relative 
to the continuum indicates that all these galaxies are 
also dominated by star formation. For the galaxies from 
Farrah et al. (2003) and Armus et al. (2007), we deter- 
mined 8 /im photometry from archival IRAC images. In 
our reductions, we included a correction for the extended 
source calibration. For those galaxies in Rigopoulou et 
al. (1999), we converted the tabulated continuum and 
aromatic 7.7 fiia peak fluxes into 8 (ith photometric val- 
ues by comparing with templates for galaxies in com- 
mon between Armus et al. (2007): IRAS 12112+0305, 
14348-1447, 1525+3609, and 22491-1808, plus Arp 220. 
The peak-to-peak scatter was a factor of two (in agree- 
ment with the finding of Chary & Elbaz (2001)). We 
add a variety of LIRGS (Alonso-Herrero et al. 2006b), 
again avoiding any galaxies with indications of signifi- 
cant levels of AGN power. We took IRAS measurements 
from Sanders et al. (2003), and computed L(TIR) to be 
consistent with the approach taken in that paper. We 
define L(8 /im) as being proportional to I'fu, and use 
similar definitions for all luminosities at specific wave- 
lengths. Figure 7 shows the relation between rest-frame 
L(8 /im) and L(TIR) for these z r^ galaxies. 

4.4.2. Results from the literature at z ^ 0.85 

We add a sample of infrared galaxies at z ^ 0.85 from 
Marcillac et al. (2006). This paper utilizes 15 /im data 
from ISO as well as 24 /im data from Spitzer and, for a 
subset, radio data from the VLA; the authors found that 
the galaxies had rest-frame L(8 /xm)/L(12 /im) ratios at 
the high end of the local galaxy dispersion. We make 
the assumption that L(12 /im) correlates with L(TIR) in 
that sample, as it is observed to do well at z = 0. We 
computed rest 8 /im luminosities for the members of the 
sample with 0.55 < z < 1.2 and adequate 15 /im pho- 
tometry, and corrected them to the IRAC 8 /im band by 
comparison with the IRS spectrum of IRAS 2249, a lo- 
cal star-formation-dominated ULIRG. We took L(TIR) 
from the Dale & Helou (2002) template fits (Marcillac et 
al. 2006). The results are shown in green in Figure 8. 

There is no mid-IR spectroscopy for these galaxies, but 
if they had strong mid-IR contributions from AGN we 
would expect substantial dispersion in the L(8/im) to 
L(TIR) ratio. Such a large dispersion is not seen. 

4.4.3. Results from the Literature ai z ~ 2 

We now compare the behavior of the 8/im luminosity 
vs. L(TIR) at z ^ 2 with that seen in local luminous 
galaxies. We first discuss three results from the litera- 
ture, and then show how our new measurements comple- 
ment them. 

1^ IRAS 12112+0305 and 14348-447 arc double nuclei sources. 



First, we took the stacked z « 1.9 galaxy tem- 
plate of Daddi et al. (2005) and compared the ratio 
of observed 850 /im to 24 /im flux densities, referred to 
the appropriate rest wavelengths for each of the local 
star-forming galaxies where there is sufflcient data to 
construct an accurate SED template (IRAS 1211, 1434, 
17208, 2249, and Arp 220). Excluding Arp 220, for 
which the ratio is 54, the average value for the other four 
galaxies is 14, whereas the Daddi et al. (2005) template 
indicates a value of 6.5. Thus, the template indicates 
stronger emission at rest 8 /im relative to the submm 
than is typical of local star-forming ULIRGs. 

We also compare with total infrared luminosity de- 
rived from submm measurements. Since the rest wave- 
length for the observed 850 /im point in the Daddi et al. 
template is 293 /im, we construct the average ratio of 
L(TIR)/L(293 /im) for local LIRGs and ULIRGs with 
4 X 10" < L(TIR) < 2.5 x lO^^ Lq and high quahty 
measurements at 350-450 /im (IRAS 09111, 1211, 1434, 
and 17208, Arp 220, and NGC 1614). We correct the 
flux densities for the differing effective filter bandwidths 
using the IRAS 2259 template. We correct to 293 /im as- 
suming a spectral slope of -3.3, determined from a larger 
sample of local ULIRG templates. We find 

L(TIR) = (98 ± 40) L(293 /im) 
and then plot the Daddi et al. point as a red square in 
Figure 7. We did not fit the observations of lower lumi- 
nosity galaxies, because the scatter at lower luminosities 
is very large, probably due to varying amounts of cold 
and warm dust in the infrared-emitting objects. 

Second, we make use of the stacking analysis 
carried out by Papovich et al. (2007). They re- 
port observed-frame 24 /im and 70 /im stacked photom- 
etry for galaxies that are both undetected in X-rays and 
have non-power-law IRAC SEDS, and are therefore as- 
sumed to be powered predominantly by star formation. 
We computed the rest 8 /im and 24 /im flux densities 
for this category (since they quote values for z — 2, the 
rest wavelengths are represented directly by the mea- 
surements) . We corrected the 8 /im flux density (derived 
from that observed at 24 /im) to the equivalent for the 
IRAC band, again using the SED of IRAS 2249. Finally, 
we compute L(TIR) from the relation derived from our 
local galaxy sample: 

L(TIR) = 86.5 X L(24 /im, rest)", where a = 0.9758. 

The rms scatter around this relation is only 0.14 dex. 
The resulting values are shown in Figure 7 as red dia- 
monds (we omitted the lowest flux density bin from Pa- 
povich et al. because the signal to noise ratio at 70 /im 
is too low). 

Third, we use the results of Yan et al. (2007). 
We select those galaxies that have z > 1.5, and rest- 
frame 6.2 /xm equivalent widths (measured by Sajina et 
al. 2007) greater than 0.3 /im. From Armus et al. (2007), 
this equivalent width threshold should select galaxies 
dominated by star formation. The galaxies selected 
(MIPS 289, 8493, 15928, 16144, 22530) are aU catego- 
rized as "type 1" by Yan et al. (2007), meaning their 
spectra contain strong aromatic features. The median 
rest-frame EW(7.7 /xm) is 2.5 /xm, corroborating that 
star formation dominates the mid-IR spectra. We con- 
verted the tabulated L(5.8 /im) for these galaxies into 
L(8 /im) from an average template for our local ULIRG 
sample (IRAS 2249). 
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We estimate L(TIR) for these sources using the radio- 
infrared relation. We describe our adjustment of this 
relation for use at high redshift in detail in a forthcom- 
ing publication (N. Seymour ct al., in preparation) but 
give a short summary here. Yun ct al. (2001) have quan- 
tified the relation for a large number of local galaxies, 
finding that the average ratio of flux density at 60/im 
to 1.4GHz is 110 independent of luminosity. The results 
of Yun et al. (2001) have not been K-corrected; since 
the infrared SEDs drop rapidly toward the blue, while 
the radio spectrum drops slowly, these corrections are a 
fairly steep function of redshift. They enter at a signifi- 
cant level only for the high luminosity galaxies, because 
they tend to be at higher redshift than the lower luminos- 
ity ones. Including them increases the indicated ratio to 
~ 170 at a L(TIR) = lO^^ l^ ^ second issue is that the 
radio outputs of LIRGs and ULIRGs often show signs 
of free-free absorption, which can decrease the outputs 
at 1.4GHz. To estimate this effect, we took the 19 star 
forming galaxies with adequate data from Condon et al. 
(1991) and computed an average slope for them of —0.58 
between 1.49 and 8.44 GHz. Taking the 8.44 GHz value 
to represent the un-absorbed radio output and assuming 
the intrinsic slope is —0.7 (Condon et al. 1991), we revise 
the IR/radio ratio downward to 140 to take account of 
absorption; this value is calculated by extrapolating from 
8.44 GHz to 1.4 GHz and comparing results for slopes of 
—0.58 and —0.7. The second concern is that the radio- 
infrared ratio may be a function of redshift (Kovacs et al. 
2006; Vlahakis et al. 2007). We tested this suggestion by 
determining a K-correction between 1.4 GHz and 850/im 
for an average of local LIRG and ULIRG templates (to 
be published by Alonso-Herrero et al.) and applying it 
to a large sample of high redshift sources. We find that 
the K-correction accurately reproduces the observed be- 
havior with no offset, indicating there is no significant 
evolution of the relation. We have therefore used a value 
of 140 for the intrinsic flux density ratio at rest wave- 
lengths to convert the radio measurements to TIR ones. 
We extrapolated the radio fluxes observed to a rest fre- 
quency of 1.4 GHz with a slope of —0.7. The redshift 
usefully makes the observed rest frequency significantly 
higher than 1.4 GHz, reducing the effects of free- free ab- 
sorption for the high-z galaxies. To convert to L(TIR), 
we use the relationship from our local sample, 

L(TIR)= 740 X L(60)", where a = 0.9427 
with an rms scatter of 0.08 dex. 

It is important to point out that both of the corrections 
we have derived have the effect of reducing the change we 
find in L(8^m)/L(TIR). That is, using the relation de- 
rived by Yun et al. (2001) and applying the proposed 
change at high redshift (Kovacs et al. 2006; Vlahakis 
et al. 2007) would have both resulted in a much larger 
change than reported. Therefore, we are taking a con- 
servative approach. 

We take radio flux densities for the galaxies observed 
by Yan et al. (2007) from Condon et al. (2003) as tabu- 
lated on the SSC website^°. Two of the remaining galax- 
ies are detected at low significance as indicated by in- 
spection of the grey scale images for them (MIPS 429 
and 506), so we enter rough values of 0.1 mJy for them. 
For MIPS 289, which is undetected, we assigned an up- 

^'■' http://ssc.spitzcr.caltcch.edu/fis/cxtragal/vla.html 



per limit of 0.1 mJy. We reject the source MIPS 8493 
because of fringing in the VLA image at its position. 
The other four (MIPS 289, 15928, 16144, and 22530) are 
in regions where the radio images are clean. We correct 
their observed radio measurements to the rest frame as- 
suming a radio k-correction of (1 -I- z)'''^, since a typical 
slope for a starburst radio spectrum is —0.7 (Condon et 
al. 1991). We plot the resulting L(TIR) and L(8 /zm) val- 
ues in Figure 8 as circles (filled for the detections, open 
for the upper limit.) 

4.4.4. 8 fim-to-TIR ratio for the « ~ 2 lensed galaxies 

The lensed galaxies described in this paper bridge the 
luminosity divide between the z ^ 2 extremely high- 
luminosity galaxies in the literature, and the compari- 
son LIRG and ULIRG samples at z ~ and z — 0.85. 
For the lensed sample we can determine which galaxies 
have mid-IR spectra dominated by star formation, using 
the aromatic feature strengths. For the subset of lensed 
galaxies with submm measurements and good upper lim- 
its at 70 /im, the SEDs are well enough constrained for 
a direct determination of L(TIR). Therefore, the lensed 
galaxies can cleanly test for evolution with redshift in the 
behavior of the L(8 ^m)/L(TIR) ratio. 

We use the aromatic-to-mid-IR flux ratio introduced 
in §4.2 to select galaxies whose mid-IR spectra are dom- 
inated by star formation. The ratio selects A2219a, 
Al689b, A2218a, and A1835a as all having strong aro- 
matic emission; their ratios range from 7-24 xlO~^^, 
with an average value of 16 x 10""'^'^. By comparison, 
our sample of local star-forming LIRGs and ULIRGs 
has a range in ratio of 7-32 xlO~^^, with an average 
of 21 x 10^^^. That is, the relative aromatic feature 
strength is virtually the same for the two samples, with 
no evidence that AGN have signiflcantly augmented the 
8 /im fluxes of our star-forming high-z sub-sample. The 
EW(6.2 /im) criterion we used for the Yan et al. (2007) 
sample would have returned the same sample of star- 
forming lensed galaxies, except that A2219a lacks 6.2 /im 
coverage and thus would not have been selected. 

We plot the four lensed galaxies as red triangles in 
Figure 8. 

4.5. Comparison of 8 fim-to-TIR ratio at high and low 

redshift 

The L(8 /im)/L(TIR) ratio is important for two rea- 
sons: 1) it is one of the most accessible measures of the 
physical behavior of infrared emission as a function of 
redshift; and 2) many studies rely on fluxes at observed- 
frame 15 or 24 /im to deduce L(TIR) and related proper- 
ties of galaxies at z ~ 1 or 2, and are thus relying on the 
behavior of L(8 /im)/L(TIR). Figures 7 and 8 illustrate 
the important result that this ratio is well-behaved for 
star-forming galaxies out to z ^ 2.5. In fact, this good 
behavior is a reasonably demanding test, since the var- 
ious samples at z ~ 2 are subject to different selection 
effects and have L(TIR) derived in different ways, yet 
they all follow a single trend. To recapitulate the differ- 
ent methods: the lensed galaxies are selected largely on 
ampliflcation factor, and we measured L(TIR) directly 
from template flts to rest frame mid- and far-IR and 
submm photometry. The four galaxies from Yan et al. 
(2007) were selected on the basis of colors and bright- 
ness at 24 /im, and we estimated L(TIR) from the radio- 
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infrared relation. The stacked galaxies from Papovich et 
al. (2007) were selected based on faint K-band and 24 /xm 
detections; we estimated L(TIR) from rest 24 fim mea- 
surements. The stacked galaxies from Daddi et al. (2005) 
were originally selected from BzK colors indicating they 
are at z ~ 2 and were stacked based on selections at a 
variety of wavelengths, and we used rest submillimeter 
measurements to estimate L(TIR). 

Having noted the overall uniform behavior, and the 
general similarity of this behavior independent of red- 
shift, we now evaluate the evidence from Figures 7 and 8 
for modest evolution with redshift. Specifically, there is a 
tendency for the ratio L(8/im)/L(TIR) to decrease with 
increasing luminosity for the local galaxies for L > 10^^ 
Lq, but this tendency seems reduced for the galaxies at 
z ^ 2. In fact, the high redshift galaxies appear to have 
SED behavior similar to normal local galaxies of lower 
luminosity in this regard. 

We have made two tests of the significance of this dif- 
ference. The first is based on the stacking results of 
Daddi et al. (2005) and Papovich et al. (2007). Both 
studies select the galaxies to stack on the basis of K-band 
detections, and therefore the mid- and far-infrared char- 
acteristics of the galaxies should be representative of the 
entire population of IR-active galaxies (i.e., there is no 
bias toward a specific L(8/xm)/L(TIR). To have an anal- 
ogous set of data for the local galaxies, we have aver- 
aged their L(8y^m)/L(TIR) ratios in luminosity intervals 
of lO"'^; these averages are shown as the large blue trian- 
gles in Figure 7. We have then fitted the averaged data 
points, shown as a blue line. To test whether the high 
redshift stacked results were compatible with the same 
fit, we carried out a second fit in which we included the 
three stacked points. For the galaxies above 10^^ Lq, we 
then compared x^ for the low redshift galaxies only in the 
two cases. We found that this x^ is 4.5 times larger if the 
three stacked points at z '-~^ 2 are included in the fit; that 
is, there is a strong indication that a single functional fit 
cannot be used for both the local galaxies and for those 
stacked at z ~ 2.^^ 

To determine whether a few anomalous local galax- 
ies could affect this conclusion, we repeated the proce- 
dure but using a more robust average for the luminosity- 
binned ratios. For each bin with ten or more items, we 
discarded the two highest and two lowest ratios and then 
averaged the remainder. Where there were less than ten 
items, we retained the straight average. The highest lu- 
minosity bin, which is responsible for most of the change 
in slope of the relation between the two luminosities, has 
ten objects and hence its average was recomputed. The 
results were virtually identical to those with the straight 
averages, as shown in Figure 7. 

A second test was applied to the high redshift galax- 
ies measured individually, either by Yan et al. (2007), 
or by us. In the case of the Yan et al. sample, we ex- 
pect a strong bias toward objects with a large ratio of 
L(8/xm)/L(TIR), because they selected relatively bright 
objects for spectroscopy, given the sensitivity limitations 
of the IRS. Our sample was selected on a combination of 
24/zm flux density and lensing amplification, so a bias is 
likely, although perhaps less severe. 

^^ This analysis makes the assumption that low-luminosity 
galaxies have similar behavior at low and high redshift. 



To compare with these two studies, we selected a lo- 
cal sample with relatively high 8 yum output by removing 
all the galaxies that are below the fit made to the aver- 
aged values (above). This sample is shown in Figure 8. 
It is noteworthy that the scatter in L(8/xm)/L(TIR) is 
small, with no outliers toward high L(8/im), so this is an 
appropriate comparison sample for cases where a bias is 
suspected. As before, we carried out two fits, one just 
to the local galaxies, and another to them plus the eight 
individually measured high redshift ones. We found that 
X^ for the low redshift sample with L > 10^^ L0 in- 
creased by a factor of 2.5 if we included the high redshift 
galaxies in the fit. Comparing Figure 8 with Figure 7, 
the discrepancy would be substantially larger if we had 
included all the local galaxies in our fits. This result con- 
firms that found for the stacked samples, namely that a 
single functional fit is not a good fit to the data over all 
redshifts. 

Another comparison can be made with the sample of 
infrared galaxies at z ~ 0.85 from Marcillac et al. (2006). 
The results fall in a region occupied by both low and 
high redshift fits in Figure 8, which confirms the rela- 
tively modest change in L(8y^m)/L(TIR) with redshift, 
but (not surprisingly given the small difference between 
z ~ and 2) does not definitively deviate from either fit. 
In making this comparison, we have assumed that the 
galaxies studied by Marcillac et al. (2006) have a bias 
toward bright 15/xm observed flux densities. If instead 
we treat them as an unbiased sample and average them 
in luminosity intervals, they fall along (or slightly above) 
the fit for the z ^ 2 sample in Figure 7. 

The tendency at high redshift for relatively strong 8/im 
emission could arise either through a contribution of an 
extra, AGN-powered component at 8 /im (e.g. the ex- 
tra 5.8 fxm continuum claimed by Sajina et al. 2007 and 
the niid-IR excess galaxies with hard X-ray emission dis- 
cussed by Daddi et al. 2007b), or to a change in the 
properties of the stellar-powered infrared-emitting re- 
gions. The AGN explanation runs counter to our selec- 
tion of individual galaxies with strong aromatic features 
that correspond to equivalent widths typical of purely 
star forming local galaxies. It is also counter to the pre- 
cautions taken in the stacking analyses to exclude AGN 
(e.g., Papovich et al. 2007). In addition, AGN contribu- 
tions would come in a broad range, so it would be hard 
to explain the consistent behavior of the five high red- 
shift samples in Figures 7 and 8 if AGN were causing 
the departure from the behavior from local star forming 
galaxies. Finally, the behavior we see is similar to that 
found by Zheng et al. (2007) in a stacking analysis at 
z ~ 0.7; specifically, they report that the far-IR SEDs 
appear to resemble those of local galaxies of lower lumi- 
nosity. Thus, while AGN may boost 8 nm emission in 
other sources, it is unlikely that AGN contamination is 
responsible for the behavior seen in Figures 7 and 8. 

We conclude that star-forming galaxies have a mild 
trend toward increasing L(8 /im)/L(TIR) ratio with red- 
shift that reaches up to a factor of ^ 2 at z = 2. In 
other words, L(8 /xm) does not turn over with increasing 
L(TIR) for z ~ 2 galaxies as it does at low redshift. This 
probably helps explain why some high-redshift galaxies 
have higher 6.2/im aromatic equivalent widths than are 
observed for local ULIRGs (Desai et al. 2007; Menendez- 
Delmestre et al. 2007; Sajina et al. 2007), though direct 
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comparisons cannot be made because there are not local 
galaxies with such extremely high luminosities. Because 
the lensed galaxies and stacked samples directly probe 
the lO^^-lO^^'^ Lq range at z ~ 2, we can directly com- 
pare to z ~ galaxies at similar luminosity rather than 
extrapolating. 

That aromatic luminosities do not turn over with 
L(TIR) in high rcdshift galaxies may be related to lower 
densities and optical depths of the infrared-emitting re- 
gions in the high-z galaxies. Lower optical depths would 
result in part if the metallicity were lower, resulting in a 
higher gas-to-dust ratio. Draine et al. (2007) show that 
the dust content of galaxies has a very steep dependence 
on their metallicity. The densities and optical depths of 
the star forming regions would also be lower if the high 
luminosities are associated with star formation on large 
spatial scales in the high-redshift galaxies, as opposed to 
strongly-concentrated, major-merger-associated nuclear 
star formation in local ULIRGs. Such very dense regions 
in local ULIRGs appear to account for the scattering of 
points below the z — trend line. 

Use of local templates to estimate L(TIR) for high-z 
galaxies detected at 24 ^m has indicated the existence 
of very high luminosities at z '^ 2 and beyond (e.g., 
Yan et al. 2007.) The change in L(8 /im)/L(TIR) sug- 
gested here would reduce L(TIR) for a given L(8 /im) 
by a factor of ~ 4. As shown in the figure, in our sam- 
ple there is only one galaxy selected to be dominated by 
star formation that has a L(TIR) significantly exceeding 
10^'^ Lq.^^ a luminosity of 10^"^ L0 corresponds to a 
star formation rate of r^ 2000 Mq yr~^ if we use the 
local calibration (Kennicutt 1998), and is suggestively 
close to the theoretical maximum luminosity that a star- 
burst is allowed by causality arguments — that star for- 
mation throughout a galaxy should not transpire faster 
than the dynamical timescale (Lehncrt & Heckman 1996; 
Elmcgreen 1999). (The actual maximum value depends 
on efficiency and dynamical mass.) Thus, our lower cal- 
ibration for L(8 /xm)/L(TIR) may remove the discrep- 
ancy between L(TIR) inferred for high redshift galaxies 
and plausible upper limits on what a starburst can pro- 
duce. Spatially-resolved imaging of high-redshift star- 
bursts (e.g. Nesvadba et al. 2007) can test whether they 
are in fact forming stars at the dynamical limit. 

5. SUMMARY AND IMPLICATIONS 

The 24 lira photometric band of Spitzer efficiently 
detects galaxies out to z ~ 3, and is being used as 
an extinction-robust measure of star formation rates in 
the distant universe. We present new IRS spectra and 
MIPS photometry of gravitationally-amplified galaxies, 
and use them in concert with archival Chandra images 
and published submillimeter photometry to investigate 
the nature of the very faint 24 /im sources. 

The aromatic feature flux ratios observed in our spec- 
tra agree well with those of star-forming template spec- 



tra at z '^ 0. This result implies that dust size and ioniza- 
tion distribution arc not strongly evolving with redshift. 
Thus, for z > 1 star-forming galaxies with little AGN 
contribution, the mid-IR spectra can be approximated 
by z '-~^ star-forming templates. 

At z ~ 2, the 24 /im data probe the 8 /im rest wave- 
length region. We find that the ratio of rest-frame 8 /<m 
luminosity to total infrared luminosity is approximately 
the same as for local galaxies of similar L(TIR). How- 
ever, the z '-~^ relationship does not appear to hold 
perfectly. Combining data from the literature with this 
sample, we demonstrate that z « 2 star-forming galax- 
ies have higher L(8 /im) luminosities for a given L(TIR) 
than low-redshift analogs. Functionally, the break in 
the L(8 /im) vs. L(TIR) relationship occurs at a 10- 
times-higher L(TIR) at z ~ 2 than it does at z ~ 0. 
We speculate that this is due to lower optical depths in 
the IR-emitting regions of the high-z galaxies, perhaps 
because of lower metallicity or because star formation 
is more spatially extended than in local ULIRGs (e.g., 
across a disk rather than in a compact nucleus). This re- 
sult has important implications for mid-IR-derived star 
formation rates at z ~ 2, and may indicate a maximum 
star-forming luminosity of a few xlO^^ Lq. 

The results just summarized were found for sources 
where AGN do not contribute substantially to the mid- 
IR spectra, as determined by the aromatic-to-mid-IR flux 
ratio. But 6 of the 19 mid-IR spectra in our extended 
sample do show a strong AGN. Including other AGN di- 
agnostics (X-ray detections and optical emission lines), 
8 of the 19 galaxies show evidence for nuclear activity. 
The prevalence of AGN in this sample may contradict 
the frequent assumption that the 24 /xm flux density in 
100 /xJy sources at 1 < z < 3 is dominated by star for- 
mation, with AGN a modest contaminant. 
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scope, and the Spitzer Space Telescope. This research 
has made use of NASA's Astrophysics Data System Bib- 
liographic Services and the NASA/IPAC Extragalactic 
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^^ This galaxy was selected as star— forming based on its 
EW(6.2 ^m) (Sajina et al. 2007), but its EW{7.7 ^m) is unusu- 



ally low, which may indicate a substantial AGN contribution. 



APPENDIX 

NOTES ON INDIVIDUAL SOURCES 
Source A2390a 

Johan Richard obtained a Keck spectrum for this source, shown in Figure 9. The strong Lyman a emission is 
redshifted by z = 2.858, and the line velocity spread (FWHM ^ 945 km s^^) indicates that the source is an AGN. 

Source A 226 la 

This sub- mm galaxy (Chapman et al. 2002) lacks a spectroscopic redshift; submillimeter and radio-based photometric 
estimates place it at z > 1.6 Chapman et al. (2002) and 2.5 < z < 4 (Aretxaga et al. 2003). It is the only source 
without a redshift in our sample. Our spectrum should be sensitive to the 11 /im complex for 0.9 < z < 1.8, and to the 
7.7 /im aromatic feature for 1.8 < z < 3.1. Thus, for the aromatic features to fall outside the IRS spectral coverage, 
the source would have to have very low (z < 1) or very high (z > 3.1) redshift. If the source were at z > 3.1, it must 
be ~ 1000 times more luminous than M82 to have produced the observed 24 fim flux density. It seems more probable 
that this source lies at 1 < z < 3 and lacks aromatic emission. 
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TABLE 1 
The IRS lensed sample. 





RA(J2000) 


DEC(J2000) 


sub-mm name 


ref 


PID 


z 


ref 


amp 


Our sample 


MS0451a 


04:54:07.12 


-03:00:36.4 






1 


1.95 


1 


2.5 


A851a 


09:42:54.55 


+46:58:44.3 


SMM J09429+4658 


A 


2 


2.38 


1 


1.3,2.0 


A1689a 


13:11:28.04 


-01:19:18.7 






2 


1.1537 


2 


8.5 ±2.3 


A1689b 


13:11:29.14 


-01:20:46.5 


SMM J13115-01208 


B 


2 


2.63 


2 


25-47 


A1835a 


14:01:04.96 


+02:52:24.8 


SMM J140105+025223.5 


C 


12 


2.565 


3 


3.5 + 0.5^ 


A2218a 


16:35:54.18 


+66:12:24.8 


SMM J163554.2+661225 


D 


1 


2.516 


4 


22 ±2 


A2218b 


16:35:55.16 


+66:11:50.8 


SMM J163555. 2+661150 


E 


2 


1.034 


5 


6.1 


A2218C 


16:35:59.12 


+66:11:47.4 






2 


0.97 


1 


6.7 


A2219a 


16:40:19.50 


+46:44:00.5 


SMM J16403+46440 


F 


1 


2.03 


1 


3.6 


A2261a 


17:22:21.00 


+32:07:05.5 


SMM J17223+3207 


F 


12 


l<z<3 


1 


3.3 


A2390a 


21:53:33.24 


+17:42:10.6 






12 


2.858 


6 


10.8 


A2390b 


21:53:34.43 


+17:42:21.6 






12 


0.913 


7 


3.35 


A2390C 


21:53:34.540 


+17:42:03.0 


'Source 14' 


G 


2 


0.91 


8 


10.4 


AC114a 


22:58:49.79 


-34:48:47.0 






2 


1.47 


1 


9.7 


A2667a 


23:51:40.00 


-26:04:52.0 






2 


1.034 


9 


17 


Additional sources from the literature 


Teplitz-1 


3:32:44.00 


-27:46:35.0 






8 


2.69 


10 


1.0 


Teplitz-1-BzK 


3:32:38.52 


-27:46:33.5 






8 


2.55 


10 


1.0 


Teplitz-2 


3:32:34.85 


-27:46:40.4 






8 


1.09 


10 


1.0 


A370a 


2:39:51.87 


-1:35:58.78 


SMM J02399-0136 


A 


9 


2.81 


11 


2.45 



Note. — Columns: 
1) source; 
2—3) coordinates; 

4) sub-mm name, if source has a sub-mm detection in the literature; 

5) reference for detected sub-mm counterpart: "A" is Smail et al. (2002); "B" is Knudsen (2004) and Knudsen (2007); "C" 
is Smail et al. (2000); Ivison et al. (2000); "D" is Kneib ct al. (2004); "E" is Knudsen et al. (2006); "F" is Chapman ct al. 
(2002); "G" is Cowie et al. (2002). 

6) Spitzer program in which IRS spectra were obtained: ''!" signifies PID 82; "2'' signifies PID 30775; "8" signifies PID 
252; "9" signifies PID 3241; 

7) redshift; 

8) redshift reference: "1" means redshift was determined from our IRS spectra; ''2" is based on Keck spectra (J. Richard 
et al. in prep.); "3" is from Prayer et al. (1999); "4" is from Kneib et al. (2004); "5" is from Ebbels et al. (1998); "6" is 
from Keck spectra presented in the Appendix; "7" is from Pello et al. (1991); "8" is from Frye & Broadhurst (1998); "9" 
is from Sand et al. (2005); "10" is from Teplitz ct al. (2007); "11" is from Prayer et al. (1998) and also Lutz et al. (2005). 

9) amplification estimated from the lensing models. Values from the literature are footnoted; un-footnotcd values are new; 

^ The amplification for this source has been the subject of several papers; Downes & Solomon (2003) and Motohara et al. 
(2005) argue that an intervening cluster galaxy boosts the amplification to ~ 25. Smail et al. (2005) and Smith et al. (2005) 
find that the intervening galaxy is a dwarf with small velocity dispersion and thus contributes only modestly to the total 
amplification, which they estimate as 3.5 it 0.5. We note the controversy and adopt the Smail et al. (2005) amplification. 



IRS spectra of lensed galaxies 



13 







TABLE 2 










Exposure 


TIMES. 






Cluster 


t(24 ^m) 


t(3.6 fim) 


t(IRS, 


LL) 


t(Chandra) 


MS0451 


3.68 


2.40 


3.66'' 




53 


A851 


2.78 


2.40 


6.34=- 







A1689 


2.78 


2.40 






40 


a 






7.31=^ 






b 






14.63*= 






A1835 


2.77 


6.00 


3.66*^; 


3.66=^ 


29 


A2218 


2.78 


2.40 






58 


a 






3.66'' 






b 






3.66''- 






c 






7.31^ 






A2219 


2.78 


2.40 


3.66'' 




41 


A2261 


2.77 


2.40 


5.85''; 


7.31'^ 


33 


A2390 


2.77 


2.40 






109 


a 






3.66''; 


5.61'^ 




b 






3.66''; 


5.61'' 




c 






7.31=^ 






AC114 


2.77 


2.40 


7.31=' 




70 


A2667 


2.77 


2.40 


1.95=^ 




9.2 



Note. — Exposure times, in kiloseconds, for the following bands: 
Spitzer/MIPS 24 fim; Spitzer/IRAC 3.6 pm; and Spitzer/IRS long- 
low grating. The 24 /j.m exposure time is the median (per-pixel) 
exposure time within the central 4 by 4 box. 
^ Both LLl and LL2 were given this exposure time. 
^ Only order LLl. 



TABLE 3 

Long- WAVELENGTH PHOTOMETRY AND TEMPLATE FITS. 





f(3.5 fini) 


f(4.5 fini) 


f(5.7 ^m) 


f(8.0 ^m) 


f(24 ^m) 


(PAH/tot) 


f(70 fj.m) 


L(TIR) 


template 


MS0451a 


32 ±2 


49 ±2 


52 ±3 


36 ±3 


1.32 


6.3E-11 


5.7 ±1.6 


< 2.2 X 10" 


CE ^ 


A851a 


29 ±3 


44 ±3 


58 ±3 


45 ±4 


0.63 


8.4E-11 




2.7 (2.3-3.6) xlO''^'' 


DH 


A1689a 


36 ±2 


35 ±2 


28 ±2 


66 ±4 


0.58 










A1689b 










0.32 


2.4E-10 




1.6 (0.43-2.4) xlO" 


DH 


A1835a 


88 ±2 


105 ±2 


136 ±3 


124 ±5 


0.99 


1.4E-10 


< 13 


6.8 (4.2-11) Xl0i2 


DH 


A2218a 


67 ±2 


70 ±2 


92 ±6 


105 ± 10 


1.16 


2.0E-10 


<7 


7.6 (5.7-9.5) xlO" 


DH 


A2218b 


182 ±4 


147 ±3 


107 ±3 


107 ±5 


1.67 




7.4 ±1.5 


4.8 (2.8-5.3) xlO" 


DH 


A2218C 


99 ±3 


79 ±3 


67 ±2 


523 ±11 


0.50 




10 ±1 


1.3 (0.88-2.0) xlO" 


Mrk 231 


A2219a 










0.82 


7.1E-11 




2.9 (2.7-3.4) xl0^2 


DH 


A2261a 


142 ±4 


138 ±4 


117 ± 3 


156 ±6 


0.58 




<3.2 


C 




A2390a 


39 ±3 


59 ±3 


90 ±3 


196 ±7 


0.83 


5.8E-11 


< 7.6 






A2390b 


462 ±8 


334 ±6 


258 ±5 


199 ±7 


0.88 




5.3 ±1.6 


1.7 (1.1-2.5) xlO^^ 


Mrk 231 


A2390C 










0.60 




<8 


8.2 (5.5-11) xlO"' 


DH 


AC114a 


91 ±3 


97 ±3 


78 ±3 


76 ±4 


0.41 


1.4E-10 








A2667a 


466 ±6 


349 ±5 


269 ±4 


320 ±8 


1.52 


















Literature sampl 










Teplitz-1 










0.13 










Teplitz-1-BzK 










0.20 










Teplitz-2 










0.13 










A370a 










1.36 











Note. — Columns: 1) Source name. 2-5) IRAC photometry, in /ijy. Formal errors from the photometry are quoted, 
though field crowding probably and aperture corrections probably introduce additional 5-10% errors. The photometry for 
A2218a is superior since it is multiply imaged (with different crowding). Crowding is too severe to quote photometry for 
A1689b, A2219a, and A2390c. 6) Observed 24 /xm flux density, in mjy, from DAOPHOT PSF-fitting. Errorbars from 
DAOPHOT are overly optimistic and not quoted; we use errorbars of 0.1 mJy for template fitting. 7) Aromatic-to-mid-IR 
fiux ratio, as defined in §4.2, for the z > 1.4 sample. 8) Observed 70 /im flux density, in mJy, from aperture photometry. 
Errorbars are dominated by the sky. 9) Total infrared luminosity, our best-fit and the range of acceptable templates, in Lq. 
10) Best-fitting template: Chary & Elbaz (2001) or Dale & Helou (2002). 
^ Mrk 231 template may also fit. 
'' Using amplification factor of 1.3. 
'^ L(TIR) cannot be usefully constrained because the redshift is unknown. 
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TABLE 4 
Measured Aromatic fluxes. 



6.2 a 7.7 a 8.3 a 8.6 a 11.3 cr 12.0 a 12.6 cr 13.6 cr 14.2 a 16.4 cr 17 



MS0451a 






28 


6 


4.8 


2 






4.3 


2 


A851a 


2.8 


0.6 


15.6 


0.8 


2.4 


1. 


2.7 


0.8 






A1689a 






12. 


5. 














A1689b 


6.5 


1.1 


21 


3 


3.8 


1. 


4.0 


2 






A1835a 


13 


2 


39 


10 


7.3 


2 


7.8 


2 






A2218a 


18.3 


0.9 


66 


3 


8.2 


1. 


13.5 


1 






A2218b 






65 


3 


4.8 


1.7 


16.3 


2 


20 


3 


A2218C 






14.9 


0.9 






1.5 


0.8 


1.3 


0.4 


A2219a 






19.3 


1.3 






4.6 


1 


5.6 


2 


A2390a 


1.1 


1 


12.4 


3 






5.8 


2 






A2390b 


















14 


10 


A2390C 






35 


11 






11. 


8 


17. 


12 


AC114a 


6.5 


0.5 


23 


2 






7.1 


1 


9.2 


1 


A2667a 






102 


11 


12 


3 


22 


4 


29 


7 



5.0 



3.8 1.2 10.8 1.3 1.7 0.8 14 5 



6.1 2.3 
6.3 1.5 

4.2 2 

9.1 3 16 3 5.6 1.9 



Note. — Measured aromatic dust emission feature fluxes and uncertainties, measured with PAHFIT. Units are 10^ '^ erg s^ em^ . Columnn 
headers give the central wavelength, in /im, of each feature. For features comprised of multiple components (like 7.7 /^m), total flux and uncertainty 
were calculated using PAHFIT's "main feature" function. 



TABLE 5 

X-RAY FLUXES AND 30" UPPER LIMITS. 

src f(2-8 kcV) f{0.5-8 keV) 

erg s^^ cm^^ erg s~^ cm~^ 

MS0451a < 1.55 x lO^^^ < 1.30 x lO^^^ 
A851a 

A1689a < 7.9 x IQ-^-' < 4.5 x lO^^^ 

A1689b < 1.1 X 10"^* < 7.2 x 10"^^ 

A1835a < 5.8 x 10"^^ < 4.4 x lO^^^ 

A2218a < 5.2 X 10"^^ < 2.6 x lO"^^ 

A2218b < 2.8 x IQ-'^-' < 2.1 x lO^^^ 

A2218C < 2.4 X lO^^"' < 1.6 x lO"^^ 

A2219a < 3.4 X lO"^"' < 3.0 x 10"^^ 

A2261a 3.7x10-15 1.8x10-^5 

A2390a 5.1 x lO^i* 7.6 x lO^^* 

A2390b 4.4 x lO^^* 2.5 x lO^^* 

A2390C < 7.5 x lO^^"' < 5.2 x lO^^^ 

AC114a < 2.8 x lO^^^ < 2.9 x lO^^^ 

A2667a < 3.2 x lO^^"* < 4.1 x 10"^* 

Note. — X-ray limits are the 3(T sky back- 
ground plus any additional counts at the source 
position, as described in §3.5. 



TABLE 6 
Aromatic ratios for SMM J163554. 2+661225 

lines ratio SINGS Starburst 

7.7fim/6.2fiin 3.6 ± 0.2 3.6(1.3-4.8) 3.9 
7.7Atm/8.6/im 4.9 ± 0.4 5.7(4.7-9.0) 5.1 
6.2^m/8.6/im 1.36 ±0.1 1.5(1.2-3.0) 1.3 

Note. — Columns: (1) line wavelengths (^m); (2) mea- 
sured flux ratios; (3) median flux ratios and 10%— 90% 
range of variation from the SINGS sample (Smith et al. 
2007); (4) flux ratios of the average starburst template of 
Brandl et al. (2006). 
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TABLE 7 
[Ne III] 15.5 fMU AND [Ne II] 12.8 AtM. 



src [Ne II] a [Ne III] 



[Ne III]/[Nc II] 



MS0451 


- 








A851 










A1689a 


0.93 


0.26 


2.0 


0.43 


A1689b 










A2218a 










A2218b 


2.9 


0.6 


0.91 


0.31 


A2218C 


- 




- 




A2219 










A2390a 










A2390b 


2.2 


0.57 


0.90 


0.69 


A2390C 


0.83 


0.47 


- 




AC114a 


1.6 


0.29 






A1835a 










A2261a 










A2667a 


2.5 


0.49 


0.92 


0.51 


Tepl 










Tl-bzk 










Tep 2 










A370a 











2.15 



0.31 



0.41 



0.37 



0.76 



0.13 



0.33 



0.22 



Note. — Measured fine structure line fluxes, from PAHFIT, in units 
of 10~ erg s~ cm" . " ■ ■ ■ " means the wavelength regime was not 
covered; "— " means the line was covered but not detected. 



TABLE 8 
AGN Diagnostics. 



X 



vifcak aromatic 



rising 



opt 



MS0451a 




X 


A851a 






A1689a 




X 


A1689b 






A1835a 






A2218a 






A2218b 






A2218C 




X 


A2219a 






A2261a 


X 


X 


A2390a 


X 


X 


A2390b 


X 




A2390C 






AC114a 






A2667a 






Tep-1 


X 




Tep-1-BzK 






Tep-2 






A370a 


X 


X 



X 
X 



Note. — Columns: 1) source name 2) X-ray detection; 
3) low aromatic flux contribution, as defined in §4.2; 4) 
mid-IR spectrum rises strongly toward the red; 5) AGN 
lines in optical spectrum. 
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' MS0451a 





r A2390b 




Fig. 1. — Postage stamps of the lensed galaxies. Images are from the Hubble Space Telescope, Advanced Camera for Surveys, filter 
F850LP, with the following exceptions: the MS0451a image is HST ACS F814W; the A2667a and A851a images are HST WFPC2 F814W; 
and for A2261a and A2219a, the images are from the Steward Observatory 90 inch, F606W. Each postage stamp is 12" by 12". Circles of 
R = 6" are drawn to illustrate the FWHM beam of Spitzer at 24 /im. 
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Fig. 2. — IRS spectra. Lower x-axes show observed wavelength in ^m; upper x-axes show rest wavelength (/im) if redshift is known. 
Y-axes plot observed flux density (mjy). Spectra obtained for the same object in different programs are plotted separately, with program 
82 plotted first. Vertical dotted lines show the positions of expected aromatic and fine structure emission features. Note that artifacts are 
present at A > 34 fj.ni. 
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Fig. 2 — continued. 
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Fig. 3. — IRS spectrum oi z = 2.516 lensed source SMM J163554. 2+661225 behind Abell 2218 (solid thick line). Ovcrplotted are two 
z ~ spectral templates: the starburst galaxy NGC 2798 from Dale et al. (2006) (thin solid line); and the average of 13 starburst galaxies 
from Brandl ct al. (2006) (thin dashed line). Crosses show the wavelengths of known aromatic components from Smith et al. (2007). 



sr e > 

Is & 

is-"' 

o " c 



5J 



ID CL o 



X 



o- ^ ? 
Q M o 
< O M 



03 






g (t 

??! 
p re g. 

Bl g^ CO 

1 >^ o 



■ O 
~i O 



,to. 



to W 



3 

a 
3 



a 
3 



X 

3 

!^ 

a 
3 



3 



4 


■ 






■ 


2 

n 


( 






■ 


8 


. 








. 


6 


1— ( 


) — 1 




■ 


4 


. 















..# 




^ 




2 








n 


. 


. 


. 



20 40 60 80 100 

-15 -1 -2 

observed f(7.7um) (10 erg s cm ) 




X 

3 

!^ 

a 
3 

06 

a 
3 



20 40 60 80 100 



-15 -1 -2 

observed f(7.7um) (10 erg s cm ) 



10 

o 





14 





12 










a 








X 


10 


3 




!^ 




a 


8 


3 




(N 




vO 


6 


g 




3 


4 



2 ■ 



i-o-i 



20 



tS 15 ■ 



10 ■ 



' i> 






20 40 60 80 100 

-15 -1 -2 

observed f(7.7um) (10 erg s cm ) 



oq 

^ 



^ 



20 40 60 80 

-15 -1 -2 

observed f(7.7um) (10 erg s cm ) 



100 



SJ 



3: 



IT C 
^ 2. 



IRS spectra of lensed galaxies 



21 



■13 



(N 






fcJO 

u 
D 

X 

< 

Oh 



o 

fcJO 

o 



14 



15 



45 



< 44 



i 
t^ 



^ 43 

(DC 

o 



42 



42 








■'5 -14 ., .2 

log obs 0.5 - 8 keV flux (erg s cm ) 



13 




43 



44 



log rest-frame L( 10-30 keV) 



Fig. 5. — Comparison of aromatic and X-ray fluxes and luminosities. SMM J163554. 2+661225 is represented by filled symbols. 
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Fig. 6. — Photometry and L(TIR) fits. Photometry and IRS spectra (thin lines) are plotted, with flux densities and wavelengths in the 
observed frame. Overplotted are best-fit templates: Dale & Helou (2002) (thick solid line), Chary & Elbaz (2001) (thick dashed line), and 
Mrk 231 (dotted line). IRAC and HST photometry are plotted for SMM J163554.2-|-661225to demonstrate that its SED is dominated by 
star formation. 
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Fig. 6. — Photometry and L(TIR) fits, continued. 
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Fig. 7. — Comparison of L(8/jm) and L(TIR) locally and for stacked measurements at z ~ 2. The individual local galaxies are indicated 
as small blue squares, while the values averaged over L(TIR,) intervals of 0.5 dex are shown as triangles. The sizes of triangles indicate the 
uncertanties in the average values for the bins, estimated from the scatter in the values. A fit to the averages is shown as a blue line. The 
black diamonds and line show a similar fit after rejecting high and low outliers in the luminosity bins (see text). The stacked results at 
2 ~ 2 are shown as a red square (Daddi et al. 2005) and two red diamonds (Papovich et al. 2007). The errors in these results have been 
consolidated into errors in L(TIR), since the conversion to this parameter is the dominant uncertainty. The error bars show uncertainties 
by a factor of 1.5 for Daddi et al. (2005) and 1.4 for Papovich et al. (2007). The red line is a fit to the 2 ~ 2 stacked points, constrained 
to agree with the local fit at 10^^ Lq. 



E 




Fig. 8. — Comparison of L(8/im) and L(TIR) locally and for galaxies at z ~ 2 measured individually. The individual local galaxies with 
above average L(8/im)/L(TIR) are indicated as blue squares. The blue line is a fit to them. The high-z galaxies from Yan et al. (2007) are 
shown as red circles, filled for those detected at 1.4GHz and open for the undetected example. The filled triangles are the lensed galaxies 
from this work. The errors have been consolidated into L(TIR). They are shown as a factor of 1.4 for the Yan et al. (2007) galaxies, 
equivalent to the ratio of flux densities at 60^m to 1.4GHz varying from 100 to 200. For the lensed galaxies, we show errors by a factor of 
1.5, based on the range of luminosities indicated by our template fits. The red line is a fit to the z ~ 2 points, constrained to agree with 
the local fit at 10^^ Lq. The green points are individual galaxies at z ~ 0.85, from Marcillac et al. (2006) and the green line is a fit to 
those points constrained to agree with the local fit at 10^^ Lq. 
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9. — Keck/LRIS spectrum of A2390a. The strong feature at A = 4690 A is presumably Lyman a at ^ = 2.858. 



